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We report the achievement of stimulated Raman adiabatic passage (STIRAP) in the microwave 
frequency range between internal states of a Bose-Einstein condensate (BEC) magnetically trapped 
in the vicinity of an atom chip. The STIRAP protocol used in this experiment is robust to external 
perturbations as it is an adiabatic transfer, and power-efficient as it involves only resonant (or 
quasi-resonant) processes. Taking into account the effect of losses and collisions in a non-linear 
Bloch equations model, we show that the maximum transfer efficiency is obtained for non-zero 
values of the one- and two-photon detunings, which is confirmed quantitatively by our experimental 
measurements. 


I. INTRODUCTION 

Since its first demonstration as an efficient tool for ex¬ 
citing vibrational states in molecular beams , stimu¬ 
lated Raman adiabatic passage (STIRAP) [SHS] has been 
applied to various situations involving multi-level sys¬ 
tems [6], such as quantum information processing I3E], 
atomic physics [IMS], cold molecules [TSHT^ and de¬ 
terministic single-photon source [50]. As regards theo¬ 
retical studies, an extensive literature exists as well [ni- 
126] . An important advantage of this technique is that it 
permits quasi-resonant multi-photon transitions without 
populating unstable intermediate states. Moreover, as an 
adiabatic transfer method, it is relatively insensitive to 
fluctuations in the experimental parameters [4]. 

In this paper, we demonstrate the use of STIRAP 
with two microwave frequencies to perform population 
transfer between hyperfine ground states of a ®^Rb Bose- 
Einstein condensate (BEC) trapped in the vicinity of an 
atom chip. Our STIRAP starts from the highest en¬ 
ergy level in an upside-down A-system m- Because of 
collision losses within the BEC in the target state, the 
maximum transfer efficiency is obtained for non-resonant 
driving fields. This in quantitatively confirmed by a the¬ 
oretical model based on non-linear Bloch equation [28] . 

More precisely, we prepare the BEC in the 
\F = 2,mp = 2) hyperfine level of the 5 ^S'i /2 
ground state and transfer it coherently into 
\F = 2,mp = l), from which the two-photon tran¬ 
sition \F = 2 ,mF = i) <—|E = 1 , = — 1 ), which 

has very good coherence properties |5M55] . can be 
addressed. Preparing the BEC in 1^ = 2, mF = 2), 
whose magnetic moment is twice as big as other trap¬ 
pable hyperfine levels of 5 ^ 5 'i/ 2 , has the advantage of 
bigger magnetic forces for equivalent electrical power 
dissipation, allowing for example to capture atoms 
from a magneto-optical trap farther away from the 
current-carrying wires |33[ 134] . or to increase the trap 
confinement during evaporative cooling |35j . 

In order to benefit from these advantages, a reli¬ 
able protocol for coherent population transfer between 
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FIG. f. Hyperfine ground states of ®^Rb involved 

in the STIRAP process, forming a upside-down A-system. 
We transfer population from state IE = 2,mF = 2) to state 
\F = 2,mF = 1) via the unpopulated intermediate level 
\F = l,mF = 1). We write : \F = 2,mF = 2) = |ei), 
jE = 2,mF = l) = |e 2 ) and |E=l,mF = l) = \g). The 
one- and two-photon detunings are respectively defined as 
A = (5i -I- ^ 2)72 and 5 = Si — 82 . 

|F = 2,tof = 2) and |F = 2,mF = 1) is required. The 
most straightforward transfer method would be to use 
a single-photon transition in the radio-frequency range. 
However, with the typical magnetic field values used in 
our experiment (less than a few tens of Gauss), the Zee- 
man effect is nearly linear |36| . and single photon tran¬ 
sitions would spread the atoms in all (equally spaced) 
Zeeman sub-levels of the F = 2 manifold. A second pos¬ 
sibility is to use a two-photon transition in the microwave 
range, with jE = 1, wf = 1) as an intermediate level. For 
such a transition, the two-photon detuning (S = 61—62 in 
figure[^ could be arbitrarily low to maximize the transfer 
efficiency, but the one-photon detuning A = (5i -|- 62)/ ‘2 
has to remain much bigger than Hi and CI 2 (which are the 
maximum available Rabi frequencies associated with the 
two microwave fields drawn in figure in order to keep 
the population of the anti-trapping state \F = 1, = 1) 

at a reasonably low level m- When the latter condition 
is fulfilled, the two-photon transition is equivalent to a 
single-photon transition with the effective Rabi frequency 
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rieff = ^^if^ 2 /( 2 |A|) [37]. It can thus be used to imple¬ 
ment a population transfer either with a tt pulse (though 
it would imply more stringent constraints on the stabil¬ 
ity of experimental parameters) |38j or with an adiabatic 
passage across resonance (which then has to be much 
slower than the time scale l/rieff)- In comparison, the 
STIRAP protocol which we consider in this paper allows 
us to perform the same adiabatic population transfer on 
a shorter time scale, on the order 1 /y/flf -|- 

This article is organized as follows. In section [ll| we 
briefly summarise the basic principles of STIRAP, and 
describe how the usual theoretical model is adapted to 
account for losses and collisions which play a major role in 
our experiment. Section |llT| is devoted to the description 
of our experimental protocol and results. In section |TV] 
we report numerical simulations of the STIRAP process, 
in good quantitative agreement with our measurements. 
Finally, we physically discuss in section jV] the need for 
non-zero one- and two-photon detunings for maximum 
transfer efficiency. 


II. THEORETICAL FRAMEWORK 


A. Basics of the STIRAP protocol 


Let us consider the three-level upside-down A system 
of figure]^ with all the atoms initially in the highest en¬ 
ergy level |ei). The two levels |ei_ 2 ) are coupled to Ig) 
by two microwave fields with time-dependent Rabi fre¬ 
quency Oi_ 2 (t) and constant frequency detuned by 
< 5 i .2 = <-^ 1.2 — Wo 1 , 2 , where ?iwoi ,2 is the energy differ¬ 
ence between |ei, 2 ) and \g). Under the rotating wave 
approximation and after an appropriate unitary trans¬ 
formation, the Hamiltonian of the system can be written 
as [aiMlIlOl : 

iL = - {—2i5i |ei) (ei| — 2 S 2 |e 2 ) ( 62 ] 

+ [|ei) {g\ + h.c.j -I- O 2 [|e 2 ) {g\ -|- h.c.]} . (1) 

In the case where the two-photon detuning J = i5i — 
62 is zero, one of the eigenstates of H takes the form 

dlMlIiQl : 


— /o 2 ^ 02 ' ~ 

vUT + ru 


Hi 




: 




which, remarkably, does not involve state \g). In order 
to transfer the population from |ei) to |e 2 ) following the 
usual STIRAP protocol, the microwave field H 2 is turned 
on first (with Hi = 0 ), such that |<i)o) = |ei) at the begin¬ 
ning of the sequence. Then, H 2 is gradually ramped down 
to zero while Hi is ramped up to its maximum value, such 
that |$o) becomes proportional to |e 2 ) at the end of the 
sequence. If the variations in the Rabi frequencies Hi 
and H 2 are sufficiently slow, the atoms will adiabatically 
follow the change in |$o), ^md thus be transferred from 
|ei) to |e 2 ) without populating the intermediate state 


\g). More precisely, the adiabaticity condition reads |3] : 
02 < \Eo-where e = (HiH 2 -HiH 2 )/(ff?+f 22 ), 
Eq is the energy of state |<I>o) and E± are the two other 
eigenvalues of H. For example, in the particular case 
where Hi = Hq cos[7rt/(2T)] and H 2 = Hq sin[7rt/(2r)] 
with 0 < t < T (see the following sections for a possi¬ 
ble practical implementation), the adiabaticity condition 
takes the following simple form : 

r»l/(^H2-f A2-|A|), (3) 

where r is the duration of the STIRAP process. In this 
simple model, the maximum efficiency is obtained for 
zero one- and two-photon detunings (A = <5 = 0), a con¬ 
dition that will change in the following sections when 
considering a more realistic experimental situation. 


B. Modelling losses and collisions 

Inter-particle interactions play a significant role in our 
experiment, all the more that we use BECs with rela¬ 
tively high density. These effects include both energy 
shifts and collisional losses m- In the following, we will 
take these effects (and other loss mechanisms) into ac¬ 
count using the formalism of non-linear Bloch equations. 

Following |35|, we model the effects of surroundings 
upon our three-level system by introducing an ensemble 
of additional states called a bath, and by considering the 
reduced density matrix p obtained by tracing over bath 
variables. Throughout this paper, we will use the stan¬ 
dard notation pij = {l\p\j), with l,j G { 61 , 5 , 62 }. In 
the absence of loss and collisions, the evolution of p is 
described by ifidp/dt = 

The effects of atomic losses due to the anti-trapping 
behaviour of state |g) is modelled by a constant loss rate 
F. We neglected collisional effects for atoms in state |g) 
because the population of the latter remains very low at 
all times. Atoms confined in the two other states |ei) 
and 1 62 ) are subject to collisional effects. In the mean 
field theory, these can be modeled m by adding to the 
Hamiltonian Q a term of the form : nog{pe-^ei |ei) ( 6 i| -I- 
06262 |e 2 ) ( 62 !), where ng ~ 5T0^^ cm“^ is the typical 
density of our BEG, and g is the coupling constant re¬ 
lated to the scattering length a through g = ATrfi^a/m, 
with m the atomic mass and a ~ 5.77 nm for ®^Rb ED- 
Although the typical value of ngg/h in our experiment 
(on the order of a few kHz) is not negligible with re¬ 
spect to other parameters, mean field collisions shift in 
our model only lead to very small visible effect on the 
numerical simulations of the STIRAP process that will 
be described in the following sections. 

On the other hand, inelastic collisions of the form 
|2,1) -I- |2,1) —|2,0) -I- |2,2) play a significant role in 
the dynamics of our system. Note that atoms falling 
into |2,0) are untrapped and therefore lost. For this rea¬ 
son, we shall neglect the reverse collision process. Other 
processes with different output and input states are also 
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neglected, because of conservation rules. The rate of 
the relevant collision process, which can be rewritten as 
|e 2 ) + |e 2 ) -t | 2 , 0 ) + |ei), has the form TcoiPe 2 e 2 , where 
Tcoi = 722^^0 can be obtained from the value of 722 pub¬ 
lished in the literature (typically, 722 — 10 “^® m^.s“^ 
[42] . which results in V^oi — 50 s“^ in our case). Follow¬ 
ing [55], we model this collision process by the combina¬ 
tion of a loss process on level |e 2 ) and a population flow 
from |e 2 ) to |ei), both occurring at the rate Ti.oiPe 2 e 2 - 
This procedure leads |28| to the following equations for 
the population terms of the reduced density operator : 


Peiei 2 Pgei) ~‘^^colPe2e2 i 

_ '^1 'l _l_ -^2 / N 

Paa ~ * 2 P<iia) * 2 P<i2g) 


^Pgg ^ 


Pe2e2 “ * 2 ^P^^a ^962) 2rcoZPe2e2 ’ 


(4) 


and for the coherences : 

— -A , .^2 I '^1 / '1 

Peig — l^lPeig 4 ” * 2 4 ” * 2 Paa) 


I Pei 


Oo 


Peie 2 - ^ (^1 ^ 2 ) 

rcoiPe2e2Peie2 > 


.^1 


Peie 2 4" * 2 * 2 


_ .^2 , , .f^l .J 

Pge2 — \Pgg ~ ^6262) ~ *^Peie2 ~ 102 Pge 2 

_£ _r 

2^962 ^ colPe2e2Pge2 j 


(5) 


where (5i = i5i - nogpeiei and 82 = 82 - noppe 2 e 2 in¬ 
clude mean field effects. This set of equations will be 
the starting point for the numerical simulation of our ex¬ 
perimental sequence, which will be described in the next 
sections. 



FIG. 2. Relative amplitude of the microwave fields a and 
/? used for the STIRAP protocol, as a function of time; solid 
blue line and open blue circles : a oc tuned to the 

transition between |ei) and 1^); red dashed line and filled 
red circles : P oc 122 (t) tuned to the transition between |e2) 
and Ip). Circles are the measured values at the output of 
the amplifier and solid and dashed lines are plotted after the 
model of equations (j^. 


where t € [0,t]. This results in a microwave with mag¬ 
netic field of the form B^vv oc (P^mwi + J 3 mw 2 )'n, where 
u is a unit (linear) polarization vector and : 

oc COs[2,7T(^fJJn^J -1- frf)i] , 

^mw2(P) OC By,.j^/3(t) COs[27r(/,22^ frf)P\ 5 ('^) 

corresponding to two side bands fmw + frf and fmw — 
frf with their relative amplitudes a(t) = sin[7rt/(2r)] 
and (3{t) = cos[ 7 rt/( 2 T)] changing over time. In order to 
take into account the non-zero extinction ratio e of our 
modulator (typically e = 2 % in our case), we introduce 
the following phenomenological expressions for a and /? : 


a{t) oc 


(1 — e) sin^ 



+ e ) 


III. EXPERIMENT 


/3{t) oc 


(1 — e) cos^ 



( 8 ) 


A. Protocol 

We start using a setup similar to [33], with a ®'^Rb 
BEG of ( 6,8 ± 0, 6 ) • 10^ atoms, which is transferred after 
the evaporative cooling stage into a harmonic trap with 
trapping frequencies 120 Hz and 210 Hz in the horizontal 
plane, and 228 Hz in the vertical plane. The STIRAP is 
performed in the latter trap. 

The two microwave fields needed for the experimen¬ 
tal protocol are generated by I/Q modulation. For 
this purpose, we use a microwave carrier wave of the 
form C oc cos(2 tt, and two radio-frequency “in- 
phase” and “quadrature” signal waves of the form S± oc 
cos [2TTfrft + (^±(t)], with 



In practice, r is the duration of the STIRAP sequence, 
wi and UJ 2 are respectively equal to 27r(/mui + frf) and 
— frf), and Hi and 172 are proportional to a and 
p. The overall microwave signal is amplified to 40 dBm 
before being radiated by a microwave horn in the direc¬ 
tion of the atomic cloud. The modelled variations of a 
and P are shown by solid and dashed lines in figure 
The corresponding experimental power levels, measured 
at the output of the microwave amplifier, are plotted as 
circles in the same figure, showing a relatively good agree¬ 
ment. 

Since wi is mostly resonant with the transition 161)0 
\g) and 022 is mostly resonant with the transition |e 2 ) O 
|g), the Rabi frequencies 171^2 can be written as hfli = 
\/S/Apb\B±\ and hQ 2 = ^3/4^b|R|||, where ps is the 
Bohr magneton, B± = Bm^i{u — u • B/|B|) is the com¬ 
ponent of Bmwi'^ in the plane perpendicular to the lo- 


p±{t) = ± 
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cal direction of the DC magnetic field B, and By = 
Bjnw2U ■ B/|B| is the component of Bmw 2 '^ along to the 
local direction of B. For a given amount of microwave 
power, the STIRAP duration r (which must be longer 
than I/Dq to satisfy the adiabaticity condition) is mini¬ 
mized for max [f2i(t)] = max [D 2 (t)], which is expected to 
occur in our experimental setup for max = max |i?|| |, 
corresponding to an angle of 45° between u and B. 
In practice we also aim for u to be almost parallel to 
the plane of the atom chip (which is compatible with 
the latter condition), as it contains several conductive 
wires (and considering the fact that in the case of a per¬ 
fectly conducting plane the magnetic field would have to 
be parallel to the latter). Experimentally, we obtained 
Di(t = r) ~ 43 kHz and D 2 (t = 0) ~ 14 kHz (see section 


ference between these two values to residual uncertainties 
on the orientation of the horn, the magnetic field at the 
bottom of the trap and the propagation of the microwave 
field from the horn to the atoms. Still, the two values are 
sufficiently close to perform an efficient STIRAP proto¬ 
col, as will be shown in the next section. 

To distinguish between atoms in states |ei) and |e 2 ) 
and measure the transfer efficiency, we fet the atoms fall 
under gravity in a magnetic field gradient after the STI¬ 
RAP sequence, which results in their spatial separation 
as shown in figure 


IV for the estimation procedure). We attribute the dif¬ 


OD 

I 

I 0.8 

0.6 
0.4 

1 0.2 

0 

FIG. 3. a) vertical optical density (OD) profile of the BEC 
after time-of-fiight in IF = 2,mF = 2) without the STIRAP 
and b) optical density of the same BEC. c) optical density pro¬ 
file of the BEC after time-of-flight in |F = 2,mF = 1) with 
a 900 /is STIRAP sequence and d) vertical optical density 
profile of the same BEC in \F = 2, m_F = 1). Blue open cir¬ 
cle stand for experimental data, and the solid red line is a 
parabolic fit. In order to discriminate between the Zeeman 
sublevels of the F = 2 manifold, during the time-of-flight 
we apply a magnetic field gradient in the vertical direction 
of the plot. This gradient causes different accelerations for 
atoms with different mp numbers resulting in different spa¬ 
tial positions. Without the STIRAP we produce a BEC of 
(6, 8 ± 0, 6) • 10® atoms and with the STIRAP sequence we 
produces a BEC of (6, 0 ± 0, 5) • 10® atoms. The transfer effi¬ 
ciency is around 87%. 
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B. Results 

In order to determine the optimal values for the mi¬ 
crowave frequencies, we first measured the modulus of 
the magnetic field at the bottom of the trap with radio¬ 
frequency spectroscopy. We found 2,45 G ± 7 mG, which 
led, based on the Breit-Rabi formula [5S], to the fol¬ 
lowing target values : frf = 856,3 ± 2,4 kHz and 
fmw = 6,838 976 GHz ± 12,3 kHz to match the one- 
and two-photon resonance conditions. We then varied 
the three experimental parameters frf, fmw and r to 
maximize the transfer efficiency. The optimal parame¬ 
ters were found to be : r = 900 /iS, frf = 860,5 kHz 
and fmw = 6,838 945 GHz. This corresponds to a trans¬ 
fer efficiency around 87%, as illustrated in figurewhich 
shows the optical density of a BEG without and with the 
STIRAP sequence. A summary of the atomic and field 
frequencies used in the optimal case is provided in ta¬ 
ble |T] Interestingly, they correspond to non-zero values of 
the one- and two-photon detunings, a point which will 
be discussed in detail in the following sections. 

From the data in figure we observe that the max¬ 
imum optical density after the population transfer is 
higher. We attribute this to an oscillation in the size 
of the BEC, as a consequence of the rapid change (i.e. 
r lower than the inverse of the trap frequencies) in the 
magnetic potential during the STIRAP process. 

For the optimal frequencies, we plot in figure (open 
circles) the experimental transfer efficiency against the 


TABLE I. Atomic transition frequencies a;oi ,2 computed us¬ 
ing the Breit Rabi formula |36] and the measured magnetic 
field value of 2,446 G (obtained from the numerical fit, see 
section IVI, and microwave frequencies experimentally 
obtained by optimizing the STIRAP efficiency. We deduce 
the values of the two-photon detuning 5/{2'k) = 11 kHz and 
the one-photon detuning A/(27r) = -24,5 kHz. 


transition : 

|ei) ^ Iff) 

|e2) Iff) 

tJo*/( 27 r) 

[GHz] : 

6,839 824 5 

6,838 114 5 

iJil {2 tv) 

[GHz] : 

6,839 805 5 

6,838 084 5 

5 i/( 27 r) 

[kHz] : 

-19 

-30 


duration of the frequency ramp. First the efficiency in¬ 
creases to near 100 % in 900 fxs, and then it exponen¬ 
tially decreases to near 0 % at a rate of 60 s“®. For 
r = 2,5 ms [43] we plot in figure (open circles) the one- 
and two-photon resonance curves, obtained by varying 
respectively fmw (figure [^b) and frf (figure [^a). Ex¬ 
perimentally we found a 45 kHz linewidth for the one- 
photon resonance curve and a 18 kHz linewidth for the 
two-photon resonance curve. 

The dispersion in the experimental points mostly cor¬ 
responds to the shot to shot fluctuations in the atom 
number before the STIRAP sequence. 
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FIG. 4. Transfer efficiency rj as a. function of the dura¬ 
tion of the STIRAP sequence for fmw = 6,838 945 GHz and 
frf = 860,5 kHz (optimal parameters). Open blue circles : 
experimental data. Solid red line : model of equations 0 
and 0 with the fitted parameters given in section [Tv] Gyan 
surface : simulation with the fitted parameters and with noise 
as discussed in section [TV] Dash-dot magenta line : model of 
equations 0 and 0 with the fitted parameters and setting 
the one- and two-photon detunings to zero. 


IV. SIMULATION 

A. Simulation protocol 

To simulate the behaviour of our STIRAP protocol, 
we numerically integrate equations and 0 , with the 
Rabi frequency ramps of equation 0. This set of equa¬ 
tions contains five independent parameters : the mag¬ 
netic field modulus at the position of the atoms Bq = |B| 
which determines the two detunings i5i and ( 52 , the max¬ 
imum Rabi frequencies for the two microwave ramps 
^li{t = r) and 172 (t = 0 ), and the loss rates Tcoi and 
r. A fit to the three experimental data sets shown in 
figures 0 0 a and 0 b was performed to estimate these 
five parameters. It led to : r7i(t = r) = 27rx43,4 kHz, 
n 2 {t = 0) = 27rxl4,4 kHz, Bq = 2,446 G, T = 635 s"! 
and Tcoi = 160 s“^. The corresponding fitting curves are 
plotted as red lines in figures 00 a and 0 b. 

The oscillations of the simulated transfer efficiency 
with T (see red solid line in figure 0 are attributed to 
the fact that 172 (t) and r7i(0) are not exactly zero as 
they ideally should be (see equations §)• This induces a 
non-zero population of the three eigenstates of H, which 
interfere and give rise to these oscillations. 

In order to account for the noise in the current sources 
and ambient magnetic field (our experiment was per¬ 
formed without any magnetic shielding), we added in 
the simulation a magnetic Gaussian white noise, with a 
standard deviation of 2,5 mG in the time domain and an 
effective frequency range between 4,5 Hz and 450 kHz. 
About one hundred independent runs of the simulation 
were performed, with random realisations of the mag¬ 
netic field noise. In figures 0 and 0 the cyan surface 
show the two-standard-deviation dispersion of the trans¬ 


fer efficiency around its mean value. This allows us to 
capture some of the dispersion of our experimental data. 


B. Discussion of the fitted parameters 


The fitted value for I7i(7 = t) is in good agreement 
with an independent measurement of 43 kHz ± 5 kHz 
obtained by direct Rabi oscillations on the transition 
Id) -O' I 5 ). The value of Bq is also in good agreement 
with the magnetic field modulus 2,45 G ± 7 mG mea¬ 
sured by radio-frequency spectroscopy. The values of 
I7i(7 = t) and 172 (t = 0 ) are not equal for the reasons 
already discussed in section ME 

The fitted value of Tcoz is of the same order of mag¬ 


nitude as the estimated value of section |HB| For T, it 
is not straightforward to give a theoretical estimate as 
it results from multiple and complex loss mechanisms. 
However, we observe that the order of magnitude of the 
fitted value (635 s“^) is consistent with the inverse of 
tout (1-4 X 10^ s“^), defined as the time taken for an 
atom initially at rest to fall under gravity by a distance 
az = 2,3 /rm, which is half the vertical dimension of the 
BEG (obtained using the Thomas-Fermi approximation). 
Note that the repulsive magnetic force has in our case 
the same order of magnitude as the gravity force, hence 
it does not change the order of magnitude for tout- 
To conclude, our numerical model including losses and 
collisions shows a very good agreement with the exper¬ 
imental data, with realistic values of the fitted parame¬ 
ters. 


C. Numerical simulation at one- and two-photon 
resonance 

In order to emphasize the importance of having non¬ 
zero values of both the one-photon (A) and two-photon 
((5) detunings, we show in figure 0 (dash-dot magenta 
line) a numerical simulation of the transfer efficiency ver¬ 
sus T using the fitted parameters discussed above, with 
the additional condition (5 = A = 0. As expected, the 
overall efficiency of the STIRAP process is lower in this 
case. In the next sections, we will give some physical in¬ 
sights on the reason why non-zero detunings improve the 
situation in this particular case. 


V. PHYSICAL DISCUSSION: WHY DO WE 
NEED NON-ZERO ONE- AND TWO-PHOTON 
DETUNINGS? 


A. Effect of collision losses 

To make the physical effect of collision losses more 
obvious, we will use in the following a much simpler 
theoretical model than non-linear Bloch equations. For 
this purpose, we consider only a pure state |^/’(t)), and 
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fmw [kHz] - 6, 838 945 GHz 


FIG. 5. a) : variation of the transfer efficiency r] as a. function of f^f (which is related to the two-photon detuning 5 by 
5 = 4nfrf + L 0 o 2 — cuoi), for r = 2,5 ms and fmw = 6,838 945 GHz. b) : variation of the transfer efficiency rj as a fnnction of 
fmw (which is related to the one-photon detuning A by A = 2nfmw — (<^oi + a>02)/2), for r = 2,5 ms and f^f = 860,5 kHz. 
Open blue circles : experimental data, solid red line : model of equations Q and Cyan surface : simulation with the htted 
parameters and with noise as discussed in section [IV] In both sets of curve the frequency offset corresponds to the experimental 
parameters of figure]^ (optimal parameters). The vertical black dashed lines are the positions of the one- and two-photon 
resonances theoretically estimated from the htted value of the magnetic Held Bo = 2,446 G. 


we write the number of atoms in state |e 2 ) as N 2 {t) — 
I (e2|'*/’(0) where N{t) is the overall atom number 

pS] . In this framework, the instantaneous loss rate for N 2 
due to the inelastic collisions in |e 2 ) is Fcozl (e 2 |^/’(t)) P, 
leading to the following equation for N{t) : 

N = -T,oi\{e2\m)\^N . (9) 


In an ideal STIRAP sequence, [tpit)) is equal at any time 
to the eigenstate |$o(i)) (given by equation ([^) of Hq, 
the latter being defined as the Hamiltonian of equation 
0 with the additional condition i5 = 0. To go further, we 
assume that the two-photon detuning S is small enough 
to be treated as a perturbation of Hq, which is the case 
if |i5| is much smaller than all the differences between the 
eigenvalues of Hq. This condition can be rewritten as : 


l<^l 

^H2 + A2-|A| 


< 1 , 


( 10 ) 


where rig = Applying the perturbation theory to 

|$o) up to t he fi rst order in the small parameter defined 
in equation (10) then leads to a corrected state >, 

with the following expression for | < > P ■ 



(corr) 

0 


2 


/ 8SAnl\ 


( 11 ) 


Inserting the latter equation in (|^ and integrating over 
time by assuming ideal STIRAP pulses of the form 
= rig cos(7rt/(2T)) and if 2 = fig sin(7rt/(2r)) yields 
the following expression for N 2 at the end of the STI¬ 
RAP sequence : 


N 2 (r) = exp 




iV(0). 


( 12 ) 


As can be seen in this equation, the effect of collision 
losses can be mitigated by choosing a non-zero value for 
both S and A. For this purpose, the product 6 A has to be 
negative, which is indeed what we found experimentally 
when optimizing the parameters for maximum transfer 
efficiency (see table [l] and related caption). Physically, 
the reduction in losses can be understood as resulting 
from the slower growth of the population in |e 2 ) during 
the STIRAP sequence, reducing the instantaneous loss 
rate Pcoil {e 2 \ip{t)) P (at the price of a more stringent adi- 
abaticity condition ([^). In practice, |i5| cannot be made 
too big because the effect described above is counterbal¬ 
anced by losses from the anti-trapped state \g ), which can 
be shown under the same hypotheses to induce additional 
losses on N 2 of the form exp[—(5^rT/(2r2g)] [55]. 


In conclusion, we have shown in this section that the 
need for non-zero one- and two-photon detunings could 
be understand as a consequence of the collision-induced 
losses in |e 2 ). It should be kept in mind however that 
equation (12) is only valid for small values of |(5|, and 
cannot quantitatively predict the optimal values of the 
detunings in our case (although it correctly predicts the 
sign of 6 A). 


B. Effect of the imbalance of the Rabi frequencies 


In this section, we will show how the need for a non¬ 
zero one-photon detuning A is accentuated by the fact 
that the maximum values of the two Rabi frequencies 
r2i(r) and H2(0) are not equal. For this purpose, we 
proceed in the same way as for the previous section, with 
the additional hypothesis |A| ^ fig. This leads to the 
following eigenenergies for H + fiA up to the first order 
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FIG. 6. a) : variations of the eigen-energies as a function of time (equations < |13[ )), in the case <5 = 10 kHz, 172(0) = 14,4 kHz and 
f7i(T) = 43,4 kHz, for A = 0; b) : same curves with A=-6 kHz; c) : transfer efficiency (full simulation with fitted parameters) 
as a function of the microwave frequency fmw and the ratio of the maximum Rabi frequencies f22(0)/17i(r). The vertical white 
dashed line indicates A = 0. The horizontal white dashed line indicates the ratio f72(0)/f7i(r) used in our experiment. 


in |(5|/r7o and |A|/17o 
h 


En — 


E+ = 


E_ = 


02 _ q2 
A X'“2 “l 


' 


— fin 


^ 2172 


(13) 


The shape of these three energy curves versus time when 
A = 0 is shown in figure |^a, where they have been 
plotted using Rabi frequency ramps of the form fli = 
17i(r) cos(7rt/(2r)) and 172 = 172(0) sin(7r7/(2T)), with 
the following values : ^ = 10 kHz, 17i(r) = 43,4 kHz 
and 172 ( 0 ) = 14,4 kHz. The parameters used to plot 
the curves slightly violate the conditions |i5|/17o ^ 
1 and |A|/17o ^ 1 to emphasize the effect on the 
shape of the curves. With 172(0) < 17i(t) as it is 
the case in our experiment, the adiabaticity condition 
(17 i(t)172 (0 )/[t17o(7 )^] <C \Eo — E±\/h) is more stringent 
around 7 = 0. As can be seen in figurej^b, a non-zero and 
negative value for A relaxes this condition by increasing 
the minimal distance between the energy curves around 
7 = 0, making the STIRAP process more efficient. 

A more quantitative study on the interplay between 
the optimal value of A and the Rabi frequency imbalance, 
based on the full simulation of our model, is shown in 
figure l^c, where the transfer efficiency is plotted as a 
function of fmw and 172(0)/17 i(t), for 5 = 27rxll kHz. It 
can be seen, as expected, that the optimal value for |A| 
increases as 172(0)/17i(r) decreases. 

Unlike the shift in the microwave frequencies due to 
collisional losses, the shift due to imbalance of Rabi fre¬ 
quencies can be cancelled if 17i(r) and 172(0) are made 
equal. In this case, the full simulation, upper row in 
figure l^c, shows that the transfer efficiency is stilf max¬ 
imum with (5A < 0 as expected. 

VI. CONCLUSION 


We have experimentally demonstrated a microwave 
STIRAP sequence to transfer a BEC from state 


|F = 2, tof = 2) to state \E = 2, mp = 1) with a transfer 
efficiency around 87% in 900 /iS. The STIRAP optimiza¬ 
tion lead to non-zero values of the one- and two-photon 
detunings, which can be understood as a result of the mit¬ 
igation of inelastic collisions which induce losses within 
\E = 2, mp = 1). The dynamics of our experimental se¬ 
quence is described to a very good extent by a set of 
non-linear Bloch equations, with realistic values of the 
parameters. 

This transfer protocol is a useful tool for atom inter¬ 
ferometry with magnetically trapped ®^Rb, as it allows 
to prepare the atoms in \F = 2, mp = 2), where they 
have maximal magnetic moment, and coherently trans¬ 
fer them to \F = 2,mp = 1). From this state one can 
address the IE = 2,mp = 1) O \F = l,mp = —1) two- 
photon transition, which can be made very robust against 
magnetically induced decoherence. More generally, this 
work shows that microwave STIRAP between hyperfine 
ground states of magnetically trapped BECs is feasible, 
and quantitatively described by non-linear Bloch equa¬ 
tion paving the way for STIRAP-based quantum infor¬ 
mation or metrology experiments integrated on a chip. 
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